Can dark matter be a Bose-Einstein condensate? 
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We consider the possibility that the dark matter, which is required to explain the dynamics of the 
neutral hydrogen clouds at large distances from the galactic center, could be in the form of a Bose- 
Einstein condensate. To study the condensate we use the non-relativistic Gross-Pitaevskii equation. 
By introducing the Madelung representation of the wave function, we formulate the dynamics of 
the system in terms of the continuity equation and of the hydrodynamic Euler equations. Hence 
dark matter can be described as a non-relativistic, Newtonian Bose-Einstein gravitational conden- 
sate gas, whose density and pressure are related by a barotropic equation of state. In the case of a 
condensate with quartic non-linearity, the equation of state is polytropic with index n = 1. In the 
framework of the Thomas-Fermi approximation the structure of the Newtonian gravitational con- 
densate is described by the Lane-Emden equation, which can be exactly solved. General relativistic 
configurations with quartic non-linearity are studied, by numerically integrating the structure equa- 
tions. The basic parameters (mass and radius) of the Bose-Einstein condensate dark matter halos 
sensitively depend on the mass of the condensed particle and of the scattering length. To test the 
validity of the model we fit the Newtonian tangential velocity equation of the model with a sample 
of rotation curves of low surface brightness and dwarf galaxies, respectively. We find a very good 
agreement between the theoretical rotation curves and the observational data for the low surface 
brightness galaxies. The deflection of photons passing through the dark matter halos is also ana- 
lyzed, and the bending angle of light is computed. The bending angle obtained for the Bose-Einstein 
condensate is larger than that predicted by standard general relativistic and dark matter models. 
The angular radii of the Einstein rings are obtained in the small angles approximation. Therefore 
the study of the light deflection by galaxies and the gravitational lensing could discriminate between 
the Bose-Einstein condensate dark matter model and other dark matter models. 

PACS numbers: 04.50.+h, 04.20.Jb, 04.20.Cv, 95.35.+d 



I. INTRODUCTION 

The existence of the dark matter in the Universe is a well established observational fact. The rotation curves of 
spiral galaxies [lj are one of the best evidences showing a modified gravitational particle dynamic at galactic level. 
In these galaxies neutral hydrogen clouds are observed at large distances from the center, much beyond the extent 
of the luminous matter. Assuming a non-relativistic Doppler effect and emission from stable circular orbits in a 
Newtonian gravitational field, the frequency shifts in the 21 cm line hydrogen emission lines allows the measurement 
of the velocity of the clouds. Since the clouds move in circular orbits with velocity vt g {r), the orbits are maintained 
by the balance between the centrifugal acceleration vf g /r and the gravitational attraction force GM(r)/r 2 of the 
total mass M(r) contained within the orbit. This allows the expression of the mass profile of the galaxy in the form 
M(r)=rv%/G. 

Observations show that the rotational velocities increase near the center of the galaxy and then remain nearly 
constant at a value of v tgoo ~ 200 km/s [lj. This leads to a mass profile M(r) — rvf goo /G. Consequently, the mass 
within a distance r from the center of the galaxy increases linearly with r, even at large distances where very little 
luminous matter can be detected. 

Several theoretical models, based on a modification of Newton's law or of general relativity, have been pro- 
posed to explain the behavior of the galactic rotation curves. A modified gravitational potential of the form 
4> = — GM [1 + a exp (— r/ro)] / (l + a)r, with a = —0.9 and ro w 30 kpc can explain flat rotational curves for 
most of the galaxies [2[- 
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In an other model, called MOND, and proposed by Milgrom [3(, the Poisson equation for the gravitational potential 
V 2 </> = AirGp is replaced by an equation of the form V [/i (x) (|V0| /a )] — AirGp, where a is a fixed constant and 
fi (x) a function satisfying the conditions /i (x) = x for x <C 1 and /i (x) — 1 for x ^> 1. The force law, giving the 
acceleration a of a test particle becomes a — ajy for a at ^ ao and a = ^onOq for <C ao, where is the usual 
Newtonian acceleration. The rotation curves of the galaxies are predicted to be flat, and they can be calculated 
once the distribution of the baryonic matter is known. Alternative theoretical models to explain the galactic rotation 
curves have been elaborated recently by Mannheim [3], Moffat and Sokolov [j| and Roberts @. The idea that dark 
matter is a result of the bulk effects in brane world cosmological models was considered in Q ■ 

A general analysis of the possibility of an alternative gravity theory explaining the dynamics of galactic systems 
without dark matter was performed by Zhytnikov and Nester [8| . From very general assumptions about the structure 
of a relativistic gravity theory (the theory is metric, invariant under general coordinates transformation, has a good 
linear approximation, it does not possess any unusual gauge freedom and it is not a higher derivative gravity), a 
general expression for the metric to order (v/c) 2 has been derived. This allows to compare the predictions of the 
theory with various experimental data: the Newtonian limit, light deflection and retardation, rotation of galaxies and 
gravitational lensing. The general conclusion of this study is that the possibility for any gravity theory to explain the 
behavior of galaxies without dark matter is rather improbable. 

Hence one of most promising ways to explain the galactic rotation curves is by postulating the existence of some 
dark (invisible) matter, distributed in a spherical halo around the galaxies. This behavior of the galactic rotation 
curves is explained by postulating the existence of some dark (invisible) matter, distributed in a spherical halo around 
the galaxies. The dark matter is assumed to be a cold, pressure-less medium. There are many possible candidates 
for dark matter, the most popular ones being the weekly interacting massive particles (WIMP) (for a recent review 
of the particle physics aspects of dark matter see Q). Their interaction cross section with normal baryonic matter, 
while extremely small, are expected to be non-zero and we may expect to detect them directly. It has also been 
suggested that the dark matter in the Universe might be composed of superheavy particles, with mass > 10 10 GeV. 
But observational results show the dark matter can be composed of superheavy particles only if these interact weakly 
with normal matter or if their mass is above 10 15 GeV [lOj. Scalar fields or other long range coherent fields coupled 
to gravity have also intensively been used to model galactic dark matter [ll[ . 

At very low temperatures, particles in a dilute Bose gas can occupy the same quantum ground state, forming a 
Bose-Einstein (BEC) condensate, which appears as a sharp peak over a broader distribution in both coordinates and 
momentum space. The possibility to obtain quantum degenerate gases by a combination of laser and evaporative 
cooling techniques has opened several new lines of research, at the border of atomic, statistical and condensed matter 
physics (for recent reviews see [i~2l. Il3j). 

An ideal system for the experimental observation of the BEC condensation is a dilute atomic Bose gas confined in 
a trap and cooled to very low temperatures. BEC were first observed in 1995 in dilute alkali gases such as vapors 
of rubidium and sodium. In these experiments, atoms were confined in magnetic traps, evaporatively cooled down 
to a fraction of a microkelvin, left to expand by switching off the magnetic trap, and subsequently imaged with 
optical methods. A sharp peak in the velocity distribution was observed below a critical temperature, indicating that 
condensation has occurred, with the alkali atoms condensed in the same ground state. Under the typical confining 
conditions of experimental settings, BEC's are inhomogeneous, and hence condensates arise as a narrow peak not 
only in the momentum space but also in the coordinate space [14j . 

If considering only two-body, mean field interactions, a dilute Bose-Einstein gas near zero temperature can be 
modeled using a cubic non-linear Schrodingcr equation with an external potential, which is known as the Gross- 
Pitaevskii equation 

The idea that the dark matter is composed of ultra-light scalar particles (m « 10 -22 eV) in a (cold) Bose-Einstein 
condensate state, was proposed initially in [l||, and further developed in [l6j . The wave properties of the dark matter 
stabilize gravitational collapse, providing halo cores and sharply suppressing small-scale linear power. The speed of 
BEC's in atomic vapors and in galactic dark matter was studied in [T3|. A cosmological model in which the boson 
dark matter gradually condensates into dark energy was considered in [l8l |. Negative pressure associated with the 
condensate yields the accelerated expansion of the Universe and the rapid collapse of the smallest scale fluctuations 
into many black holes, which become the seeds of the first galaxies. 

Scalar mediated interactions among baryons embedded in a Bose-Einstein condensate, composed of the mediating 
particles, which extend well beyond the Compton wavelength, were studied in fl9| . If the dark matter of the Universe 
is composed of such a condensate, the imprints of an interaction between the baryonic and the dark matter could be 
manifest as anomalies in the peak structure of the Cosmic Microwave Background. In a medium composed of scalar 
particles with non-zero mass, the range of Van der Waals-type scalar mediated interactions among nucleons becomes 
infinite when the medium makes a transition to a Bose-Einstein condensed phase. In an astrophysical context this 
phenomenon was explored in [20| and the effect of a scalar dark matter background on the equilibrium of degenerate 
stars was studied. 
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A relativistic version of the Gross-Pitaevskii equation was proposed in [2l| and the cosmological implications of 
a steady slow BEC process were considered. It is interesting to note that the resulting equation of state for the 
condensate is that of the Whittaker solution [z| , which has special importance in classical general relativity. Namely, 
the solution is the static non-rotating limit of the Wahlquist solution, which is probably the most important exact 
rotating perfect fluid spacetime. 

Dolgov and Smirnov 23] assumed that the Pauli exclusion principle is violated for neutrinos, and consequently, 
neutrinos obey the Bose-Einstein statistics. Neutrinos may form cosmological Bose condensates, which accounts for all 
(or a part of) the dark matter in the universe. "Wrong" statistics of neutrinos could modify big bang nucleosynthesis, 
leading to the effective number of neutrino species smaller than three. The Pauli principle violation for neutrinos 
can be tested in the two-neutrino double beta decay. In order to effectively modify Kepler's law without changing 
standard Newtonian gravity, the galactic dark matter was described by a scalar field in [24]. For cold scalar fields, 
this model corresponds to a gravitationally confined Boson-Einstein condensate, but of galactic dimensions. A light, 
neutral vector particle associated with a vector field </> M , and which appears in a modified theory of gravity, may form 
a cold fluid of Bose-Einstein condensates before the cosmological recombination with zero pressure and zero shear 
viscosity [Hj]. Vortices in axion condensates on the galactic scale have been studied in [26|. Such vortices can occur 
as a result of global rotation of the early universe. Various cosmological implications of axion condensation have been 
investigated in [27l ] . 

The dynamical equations describing the evolution of a self-gravitating fluid of cold dark matter can be written 
in the form of a Schrodinger equation coupled to a Poisson equation, describing Newtonian gravity. It has been 
shown that, in the quasi-linear regime, the Schrodinger equation can be reduced to the exactly solvable free-particle 
Schrodinger equation. This approach can be used to study gravitational instabilities or structure formation [28| . 
Unified models for dark matter and dark energy via a single fluid and models discussing the properties of the dark 
energy and cosmological horizons have been proposed and discussed in 29] . 



It is the purpose of the present paper to systematically investigate the possibility that the dark matter, which is 
required to explain the dynamics of the neutral hydrogen clouds at large distances from the galactic center, could be 
in the form of a Bose-Einstein condensate. To study the dark matter condensate we use the non-relativistic Gross- 
Pitaevskii equation. By introducing the Madelung representation of the wave function, we formulate the dynamics 
of the (quantum) system in terms of the continuity equation and of the hydrodynamic Euler equations. Hence dark 
matter can be described as a non-relativistic, Newtonian Bose-Einstein gravitational condensate gas, whose density 
and pressure are related by a barotropic equation of state. In the case of a condensate with quartic non-linearity, the 
equation of state is polytropic with index n = 1. 

In the framework of the Thomas-Fermi approximation, the structure of the Newtonian gravitational condensate is 
described by the Lane-Emden equation, which can be exactly solved. General relativistic configurations with quartic 
non-linearity are studied, by numerically integrating the structure equations. The basic parameters (mass and radius) 
of the Bose-Einstein condensate dark matter halos sensitively depend on the mass m of the condensed particle, of the 
scattering length a and of the central density of the dark matter. The values of these parameters can be constrained 
by using the known radii and masses of the galactic dark matter halos. 

To test the viability of the model and in order to apply it to realistic systems, we compute the rotation curve of 
the Bose-Einstein condensate for a sample of 12 galaxies that includes low surface brightness and dwarf galaxies with 
measured rotation curves extending in the dark matter dominated region. A best fit to the rotation curves of galaxies 
is obtained in terms of a parametric core baryonic mass distribution, which is superposed with the Bose-Einstein 
condensate mass distribution. It turns out that the Newtonian potential of the core is asymptotically decreasing, but 
the corrected rotation curve is much higher than the Newtonian one, thus offering the possibility to fit the rotation 
curves. For low surface brightness galaxies we find a very good agreement between the theoretical rotation curve with 
a normal baryonic core and the observational data. 

The deflection of photons passing through the dark matter halos is also analyzed, and the bending angle of light 
is computed. The bending angle obtained for the Bose-Einstein condensate is larger than that predicted by standard 
general relativistic and dark matter models. The angular radii of the Einstein rings are obtained in the small angles 
approximation. Therefore the study of the light deflection by galaxies and the gravitational lensing provides a powerful 
observational method to discriminate between the Bose-Einstein condensate dark matter model and other dark matter 
models. 

The present paper is organized as follows. The basic equations describing the gravitationally bounded Bose-Einstein 
condensate are derived in Section II. The case of the condensate with quartic non-linearity is considered in Section 
III. Newtonian dark matter condensate models are analyzed in Section IV. General relativistic condensate models 
are studied numerically in Section V. We compare the predictions of our model with the observed galactic rotation 
curves in Section VI. The bending of light by dark matter condensate halos is considered in Section VII. We discuss 
and conclude our results in Section VIII. 
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II. THE GROSS-PITAEVSKII EQUATION FOR GRAVITATIONALLY BOUNDED BOSE- EINSTEIN 

CONDENSATES 



In a quantum system of N interacting condensed bosons most of the bosons lie in the same single-particle quantum 
state. For a system consisting of a large number of particles, the calculation of the ground state of the system with 
the direct use of the Hamiltonian is impracticable, due to the high computational cost. However, the use of some 
approximate methods can lead to a significant simplification of the formalism. One such approach is the mean field 
description of the condensate, which is based on the idea of separating out the condensate contribution to the bosonic 
field operator. We also assume that in a medium composed of scalar particles with non-zero mass, the range of Van 
der Waals-type scalar mediated interactions among nucleons becomes infinite, when the medium makes a transition 
to a Bose-Einstein condensed phase. 

The many-body Hamiltonian describing the interacting bosons confined by an external potential V ext is given, in 
the second quantization, by 



H 



dr^ + (f) 



h 2 

2m 



* (r) + 7: / drdr 1 ^ (r) (/) V (r - f) W (r) (r") , (1) 



where if? (r) and *f? + (r) are the boson field operators that annihilate and create a particle at the position f, respectively, 
and V (r — r 1 ) is the two-body interatomic potential [12|. V ro t (r) is the potential associated to the rotation of the 
condensate, and is given by 



V rot (r) = frot (*) ■ 



(2) 



where uj is the angular velocity of the condensate and f ro t (t) a function which takes into account the possible time 
variation of the rotation potential. For a system consisting of a large number of particles, the calculation of the ground 
state of the system with the direct use of Eq. (fT|) is impracticable, due to the high computational cost. 

Therefore, the use of some approximate methods can lead to a significant simplification of the formalism. One such 
approach is the mean field description of the condensate, which is based on the idea of separating out the condensate 
contribution to the bosonic field operator. For a uniform gas in a volume V, BEC occurs in the single particle state 
= fW~, having zero momentum. The field operator can be decomposed then in the form ^? (f) = yj N/V + \f?' (f). 
By treating the operator \f?' (r) as a small perturbation, one can develop the first order theory for the excitations of 
the interacting Bose gases 

In the general case of a non-uniform and time-dependent configuration the field operator in the Heisenberg repre- 
sentation is given by 



= V(r,t) + *' {r,t), 



(3) 



where i/j(r,t), also called the condensate wave function, is the expectation value of the field operator, ip(r,t) = 
(r, It is a classical field and its absolute value fixes the number density of the condensate through p(r,t) — 

\ip (r, t) | 2 . The normalization condition is N = J p (r, t) d 3 r, where iV is the total number of particles in the condensate. 

The equation of motion for the condensate wave function is given by the Heisenberg equation corresponding to the 
many-body Hamiltonian given by Eq. (TT]), 



<f?,H 



~2m~ V + Vrot ^ + Vext (f) + ' ( ' l) V ^-r)* (r, t) 



*(r,t) 



(4) 



By replacing ^f! (r, t) with the condensate wave function ip gives the zeroth-order approximation to the Heisenberg 
equation. In the integral containing the particle-particle interaction V (f* — r) this replacement is in general a poor 
approximation for short distances. However, in a dilute and cold gas, only binary collisions at low energy are relevant 
and these collisions are characterized by a single parameter, the s-wave scattering length, independently of the details 
of the two-body potential. Therefore, one can replace V (f — f) with an effective interaction V (f — r) = AS (r 1 — r), 
where the coupling constant A is related to the scattering length a through A = Airft^a/m. Hence, we assume that in a 
medium composed of scalar particles with non-zero mass, the range of Van der Waals-type scalar mediated interactions 
among nucleons becomes infinite, when the medium makes a transition to a Bose-Einstein condensed phase. 

With the use of the effective potential the integral in the bracket of Eq. (jU gives X\ip(r, t)\ 2 , and the resulting 
equation is the Schrodinger equation with a quartic nonlinear term (l2l |. 
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However, in order to obtain a more general description of the Bose-Einstein condensate stars, we shall assume an 



arbitrary non- linear term g \ \ip (r, t)|' 



where we have denoted 



(5) 



The probability density p is normalized according to / d 3 fp = N. As was pointed out in 31], the Gross-Pitaevskii 
approximation is a long-wavelength theory widely used to describe a variety of properties of dilute Bose condensates, 
but for short-ranged repulsive interactions this theory fails in low dimensions, and some essential modifications of the 
theory are necessary. 

Therefore the generalized Gross-Pitaevskii equation describing a gravitationally trapped rotating Bose-Einstein 
condensate is given by 



, d , . 



2m 



V r , 



(r) + V ext (r)+g' (\^(r,t)f 



i> (r, t) , 



(6) 



where we denoted <?' = dg/dp. The mass of the condensed particles is denoted by m. 

As for V ex t (r), we assume that it is the gravitational potential V, Vext — V, and it satisfies the Poisson equation 



W = 4irGp m , 



(7) 



where p m — mp is the mass density inside the Bose-Einstein condensate. 

The physical properties of a Bose-Einstein condensate described by the generalized Gross-Pitaevskii equation given 
by Eq. ^ can be understood much easily by using the so-called Madelung representation of the wave function [33 |. 
which consist in writing ip in the form 



ip{f,t) = V P (r, t) exp 



(8) 



where the function S (r, t) has the dimensions of an action. 

By substituting the above expression of the wave function into Eq. ([5]) it decouples into a system of two differential 
equations for the real functions p and v, given by 



dp„ 
dt 



+ V • (p m v) = 0, 



(9) 



dv 
di 



-VP 



\ m 1 



V r , 



p m y 



where we have introduced the quantum potential 



the velocity of the quantum fluid 



h 2 v 2 Vp 

2m yfp 



Til 



respectively, and we denoted 



P 



/ Pm\ _ i (Pm\ Pm _ / Pm\ 
V 771 / ' \777/777 \ 171 ' 



From its definition it follows that the velocity field is irrotational, satisfying the condition V x v = 0. 
The quantum potential Vq has the property 



(10) 



(11) 



(12) 



(13) 



pViV Q = V, f -— pViVj lnpj = Vjog, 

where cr^ = — h 2 p\7 i \7j lnp/4m is the quantum stress tensor, which has the dimensions of a pressure 
intrinsically anisotropic quantum contribution to the equations of motion. 



(14) 
and is an 
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Therefore the equations of motion of the gravitational ideal Bose-Einstein condensate take the form of the equation 
of continuity and of the hydrodynamic Euler equations. The Bose-Einstein gravitational condensate can be described 
as a gas whose density and pressure are related by a barotropic equation of state [3(3|. The explicit form of this 
equation depends on the form of the non-linearity term g. 

For a static ideal condensate, v = 0. In this case from Eq. (flQ|) we obtain 



Vq + Vrot + V ext + g' = constant (15) 
By applying the operator V 2 to both sides of Eq. (fl5)) gives 

V 2 (V Q + V rot + 9') + V 2 V ext = 0. (16) 

In the case of a condensate with a non-linearity of the form g (p m ) = kp 2 m j2 and in the presence of a confining 
gravitational field V — Vext, it follows that the generalized potential V gen = — Vq — V ro t — kp m satisfies the Poisson 
equation, 

W ge „ = ^G Pm . (17) 

If the quantum potential can be neglected, the mass density of the condensate is described by the Helmholtz type 
equation 

III. STATIC AND SLOWLY ROTATING NEWTONIAN BOSE-EINSTEIN CONDENSATES 

When the number of particles in the gravitationally bounded Bosc-Einstcin condensate becomes large enough, 
the quantum pressure term makes a significant contribution only near the boundary of the condensate. Hence it is 
much smaller than the non-linear interaction term. Thus the quantum stress term in the equation of motion of the 
condensate can be neglected. This is the Thomas-Fermi approximation, which has been extensively used for the study 
of the Bose-Einstein condensates 1211 . A s the number of particles in the condensate becomes infinite, the Thomas- 
Fermi approximation becomes exact 33] . This approximation also corresponds to the classical limit of the theory (it 
corresponds to neglecting all terms with powers of H or as the regime of strong repulsive interactions among particles. 
From a mathematical point of view the Thomas-Fermi approximation corresponds to neglecting in the equation of 
motion all terms containing Vp and VS. 

In the case of a static Bose-Einstein condensate, all physical quantities are independent of time. Moreover, in the 
first approximation we also neglect the rotation of the condensate, taking V ro t = 0. Therefore the equations describing 
the static Bose-Einstein condensate in a gravitational field with potential V take the form 



^P(—)=- Pm v(-), (19) 
V m I \m J 



V 2 V = AirGpm. (20) 

These equations must be integrated together with the equation of state P = P (p m ), which follows from Eq. (|13p . 
and some appropriately chosen boundary conditions. By assuming that the non-linearity in the Gross-Pitaevskii 
equation is of the form 

g(p) = ap T , (21) 

where a and T are positive constants, a — constant > 0, T — constant > 0, it follows that the equation of state of the 
gravitational Bose-Einstein condensate is the polytropic equation of state, 

P(p m )=a(T-l)p^ = Kp r m , (22) 

where we denoted K — a (T — 1). 

By representing T in the form T = 1 + 1/n, where n is the polytropic index, it follows that the structure of the 
static gravitationally bounded Bose-Einstein condensate is described by the Lane-Emden equation, 

1 d ( nd6\ 

£ 2 — +0"=O, (23) 
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where 8 is a dimensionless variable defined via p 

r 1/n-l -| 1/2 

transformation r = (n + l)X/9 C m /4irG 



Pcm,O n and £ is a dimensionless coordinate introduced via the 
£■ Pcm is the central density of the condensate [34j |. 



Hence the mass and the radius of the condensate are given by 



R = 



(n+l)a 



AnGn 



1/2 



„(l-n)/2 
rem 



6, 



and 



M = 4tt 



(n + 1) a 



4ttGV 



3/2 



(3-n)/2n^2 



e?i^(ei)i, 



respectively, where £i defines the zero- pressure and zero-density surface of the condensate, (£i ) = [34 
In the standard approach to the Bose-Einstein condensates, the non-linearity term g is given by 



/ \ "0 , , ,4 U0 2 

s(p) = yW = T^' 

where Mq = ^nY?ajm [T3 |. The corresponding equation of state of the condensate is 



-P (pm) = C^oP : 



with 



^0 = 



2nh 2 a 



(24) 



(25) 



(26) 



(27) 



(28) 



Therefore the equation of state of the Bose-Einstein condensate with quartic non-linearity is a polytrope with index 
n = 1. In this case the solution of the Lane-Emden equation can be obtained in an analytical form, and the solution 
satisfying the boundary condition (0) = 1 is [34[ 



o(0 = 



sin£ 



The radius of the gravitationally bounded system is defined by the condition 6 (£i) 
the radius R of the Bose-Einstein condensate is given by 



R = 7T 



Gm 3 ' 



(29) 

0, giving £i = vr. Therefore 
(30) 



The radius of the gravitationally bounded Bose-Einstein condensate with quartic non-linearity is independent of 
the central density and the mass of the system, and depends only on the physical characteristics of the condensate. 
The total mass of the condensate is obtained as 



M = 47T d 



h 2 a 
Gm 3 



3/2 



Pcm\0' Kl)|=47r' 



\Gm 3 J 



3/2 



(31) 



where we have used \9' (£i)| = 1/tt. 

The mass-radius relation for the static condensate is given by 



M = -R 3 p c 



(32) 



which shows that the mean density of the condensate (p„ 
the condensate by the relation 



3M/4ttR 3 can be obtained from the central density of 



/ \ Q Pcm 
Pm = 3—2-. 



(33) 



The case of the slowly rotating Bose-Einstein condensates can also be straightforwardly analyzed, by taking into 
account the fact that the condensate obeys a polytropic equation of state. The study of the slowly rotating polytropes 
was performed in detail by Chandrasekhar [351 ] - 
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The Lane-Emden equation for a rotating Bose-Einstein condensate is 



i n ( ? an 



1 d 



(1 - M 2 ) 



de 

dp 



= + a, 



(34) 



where /i = cos0 and = uj 2 /2irGp cm . The radius i? w and the mass M u of the condensate in slow rotation are given 
in the first order in f2 by 



and 



= 4tt 



(n + l)a 



47rGn 



1/2 



,(l-n)/2£ 
cm Si 



(n+ l)a 



3/2 



(3-n)/2 n;C 2 
rem 



e 2 |0'(a)i 



3^0 (6) J 1/3 
ai^'(6)i . ' 



(35) 



(36) 



respectively. The values of the function ■i/'o are tabulated in [35 1. 

In the case of the Bose-Einstein condensate with quartic non-linearity with polytropic index n = 1, the Lane-Emden 
equation can be integrated exactly, giving for the radius and mass of the rotating condensate the following simple 
relations 



R u =#(1 + 30) 



1/3 



(37) 



M,, = M 



1 



i n 



With respect to a scaling of the parameters m, a and p cm of the form m — * aim, a 
radius and the mass of the condensate have the following scaling properties: 



R 



a 



-3/2 1/2 



R,M 



-9/2 3/2 



0% a s M. 



a 2 a, p c 



(38) 



a3Pcm, the 



(39) 



IV. DARK MATTER AS A BOSE-EINSTEIN CONDENSATE 

In the present Section we analyze, by using the results obtained in the previous Sections, the possibility that dark 
matter is in the form of a Bose-Einstein condensate. For simplicity we shall restrict our study to the case of a 
Bose-Einstein condensate with a quartic non-linearity. As a first result one has to point out that, if dark matter is a 
Bose-Einstein condensate, it cannot be pressure-less, as usually considered in most of the investigations of the galactic 
dark matter problem, but must obey a polytropic equation of state with index n = 1. 

The possibility that dark matter has a substantial amounts of pressure, comparable in magnitude to the energy 
density, has been discussed in [36| . Galaxy halo models, consistent with observations of flat rotation curves, are 
possible for a variety of equations of state with anisotropic pressures. However, in the case of dark matter in the form 
of the Bose-Einstein condensate the pressure distribution is isotropic. 

The density distribution of the dark matter Bose-Einstein condensate follows from Eq. (f2"9"|) and is given by 



< \ ( c ) 
PDM (r) = p y D ' M 



sin kr 
kr 



(40) 



where k = Gm 3 /h 2 a and p^ M is the central density of the condensate, Pp M = Pdm (0). The mass profile of the 
dark condensate galactic halo ttidm {t) = 47r J q pdm {t) r 2 dr is 



m DM {r) = rf^-r 



sin kr 
kr 



cos kr 



(41) 



Eq. (|4Tj) allows to represent the tangential velocity v \ g (r) — Gm dm (r) /r of a test particle moving in the dark halo 



2 / x _ ^GPDM 
v tg V) - k 2 



sin kr 
kr 



cos kr 



(42) 
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For r-iOwe have v\ g (r) — > 0. 

At the boundary of the halo, which is defined by the radius Rdm of the halo, the density of the Bose-Einstein 
condensate is negligible, pdm{Rdm) = and kr — > 7T. These conditions give the radius of the dark matter condensate 
as 



, h 2 a 

Rdm=k\1-q^- (43) 



The total mass of the condensate is given by 



M DM = m DM {Rdm) = -RdmPdm- ( 44 ) 

Near the vacuum boundary of the condensate the maximum tangential velocity of a test particle tends to a constant 
value which can be expressed as 

„2 {R) = 4 * G rf>M = ^£dm3m (45) 
rC 7T 

Eqs. I|40 p -(|45 p give a complete description of the physical properties of the galactic dark matter condensates in 
terms of three quantities, the mass m of the condensed particle, the inter-particle scattering length a, and the central 
density of the condensate. Together with Eqs. (|30p and (f3"2")) they can be used to constrain the physical properties of 
the condensate. 

As a toy-model we consider the example of a galactic dark matter halo extending to up to Rdm = 10 kpc ~ 
3.08 x 10 22 cm, with a mass of the order of Mum = 3 x 10 11 M q . The mean density of the condensate dark matter 
is (pdm) = 3Mdm/4ttRjj M ps 5.30 x 10~ 24 g/cm 3 . Then Eq. f3"5]) gives the central density of the condensate as 

Pdm = (Pdm) /3 « 1.74 x 1(T 23 g/cm 3 . 

An important parameter in the description of the Bose-Einstein condensate is the scattering length a. In terrestrial 
experiments performed with 87 Rb a has a value of the order of a ( 87 Rb) ps 5.77 x 10~ 7 cm pa 5.77 x 10 6 fm [l4| . 
However, for very light particles a may be much larger. 

The mass of the particle in the condensate can be obtained from the radius of the dark matter halo as 

2fe2 \ 1/3 

J m 2.58 x 1CT 30 [a (cm)] 1/3 [R (kpc)p 2/3 g « 6.73 x 10~ 2 [a(fm)] 1/3 [R (kpc)p 2/3 eV. (46) 

From this equation it follows that the mass of the particle forming the condensate dark matter halos is of the order 
of eV. For a ps 1 fm and R ps 10 kpc, the mass may be of the order of m ps 14 meV. For values of a of the order of 
a ps 10 6 fm, corresponding to the values of a observed in terrestrial laboratory experiments, m ps 1.44 eV. Of course 
much larger values of a may bring the mass of the condensate particle in the range of a few tenth of eV's. 

With the use of Eq. ([45]) one can estimate the tangential velocity of a test particle in the dark matter halo. By 
using the numerical values from the previous estimations we obtain U{ S (-Rdm) ~ 365 km/s, a value which is consistent 
with the observations of the galactic rotation curves. From the definition of the mean density it follows that the total 
mass of the condensate can be obtained from the tangential velocity as Mum ~ v\ '„(Rdm)Rdm / 'G. 

The density profile of the Newtonian dark matter condensate and its mass profile are represented in Fig. [T] 

The tangential velocity of a test particle in the gravitationally bounded galactic condensate is shown, for different 
values of the central density, in Fig. [51 



V. THE GENERAL RELATIVISTIC BOSE-EINSTEIN CONDENSATE GALACTIC HALOS 

In the previous Sections we have considered the gravitationally bounded Bose-Einstein condensate dark matter in 
the framework of Newtonian gravity. General relativistic effects may change the physical properties of the dark matter 
condensate halos in both a qualitative and quantitative way. These effects may become particularly important since 
as a condensate dark matter cannot be described as a pressure-less fluid. Therefore the study of the general relativistic 
Bose-Einstein condensates offers a better understanding of their physical properties, and opens the possibility of direct 
observational testing of their properties. 

For a static, spherically symmetric distribution of matter in the galactic halo the line element is 



ds 2 = _ e u{r) c 2 dj 2 + e X(r) dr 2 + r 2 ^2 + ^2 ^2) 



(47) 
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FIG. 1: Density profile (in units of 10" 25 g/cm 3 ) (left fi gure) and mass distribution (in units of 1O 1O M0) (right figure) of the 
Newtonian Bose-Einstein condensate dark matter for different values of the central density: p^ M — 2 x 1CP 25 g/cm 3 (solid 
curve), pp M = 5 x 10~ 25 g/cm 3 (dotted curve), p ( p M = 8x 10~ 25 g/cm 3 (short dashed curve), pp M = xl0~ 24 g/cm 3 (long 
dashed curve) and p^ M = 2 x 10~ 24 g/cm 3 (ultra-long dashed curve). 
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r ( kpc ) 

FIG. 2: Tangential velocity of a test particle in a galactic Newtonian Bose-Einstein condensate for different values of the central 
density: p^ M = 2 x 1CT 25 g/cm 3 (solid curve), p^ M = 5 x 1CT 25 g/cm 3 (dotted curve), p$ M = 8 x 1(T 25 g/cm 3 (short dashed 
curve), p^ M = xlCF 24 g/cm 3 (long dashed curve) and Pp M = 2 x 1CF 24 g/cm 3 (ultra-long dashed curve). 



The galactic rotation curves provide the most direct method of analyzing the gravitational field inside a spiral 
galaxy. The rotation curves have been determined for a great number of spiral galaxies. They are obtained by 
measuring the frequency shifts z of the light emitted from stars and from the 21-cm radiation emission from the 
neutral gas clouds [l[. As shown in Appendix [X[ the tangential velocity of test particles in stable circular orbits is 
given by :37] 



(48) 



Thus, the rotational velocity of the test body is determined by the metric coefficient exp (y) only. 

A general relativistic static dark matter distribution with energy density pdm{t) and pressure Pdm{t) can be 
described by the mass continuity equation and by the Tolman-Oppcnhcimcr-Volkoff (TOV) equation, which are given 
by® 



(49) 
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dP DM (r) _ (G/c 2 ) (p DM c 2 + P DM ) (AnP DM r 3 /c 2 + m DM ) 

, (5U) 



dr 



i _ 2Gm DM (r) 

1 „9„ 



du__ 2P' DM (r) _ 2 (G/c 2 ) (ATrP DM r 3 /c 2 + m DM ) 

dr ~ pdmc 2 + Pdm ~ r ~ r " 1 ' ' ' ' 



2Gm DM (r) 



where mpM^) is the mass inside radius r. The system of equations (|4"5 ]) -(|5T) |) must be closed by choosing the equation 
of state for the thermodynamic pressure of the dark matter, 

Pdm = Pdm (pdm) ■ (52) 

At the center of the dark matter distribution the mass must satisfy the boundary condition 

m DM (0) = 0. (53) 

For the thermodynamic pressure of the dark matter Pdm we assume that it vanishes on the surface of the dark 
halo, Pdm(R) = 0. The exterior of the Bose-Einstein condensate halo is characterized by the Schwarzschild metric, 
describing the vacuum outside the galaxy, and given by [38| 

(e v r l = (e- x ) ext = 1 - 2GMdm , r > R, (54) 

where Mum = "Jj)m(-R) is the total mass of the dark halo. The interior solution must match with the exterior 
solution on the vacuum boundary of the star. 

As for the equation of state of the dark matter we adopt the equation of state corresponding to the Bose-Einstein 
condensate with quartic nonlincarity, given by 

Pdm (pdm) = U p 2 DM , (55) 

with Uo = 2TTh 2 a/m 3 . 

The structure equations for the Bose-Einstein condensate stars can be written in a dimensionless form, by intro- 
ducing a set of dimensionless variables 77, Mo, 6 and S, respectively, and defined as 

r = r*r), m DM = m*M , p DM = Pom^i p dm = p { D M c 2 \9 2 , (56) 
where p DM is the central (r = 0) value of the energy density of the Bose-Einstein condensate. By taking 

(57) 



^Gp^ M 



and 



^(r*) a p£k= ; (58) 



M 



respectively, and by denoting 



TJ n (c) 

A = —LEmAL ; (59 ) 



the structure equations of the Bose-Einstein condensate galactic dark matter halos with quartic non-linearity can be 
written as 

dM n 2 



dij 

d9 _ (1 + X o 0) (A 6>V + M ) 



dV 2 \ oV 2 (l - 2fa 



= 0»r, (60) 

(61) 
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dv 
drj 



4A d9 2 (A o 0V + M ) 



1 + A 6> d?7 



1 - 



2i\/< 
»7 



(62) 



In the new dimensionless variable the boundary conditions at the center and surface of the condensate halos are 
given by 



M o (0) = 0, 6(0) = 1, d( vs ) = 0, 



(63) 



where 775 is the value of the dimensionless radial coordinate at the vacuum boundary of the galaxy. 

The general relativistic density and mass profiles of the Bose-Einstein condensate dark halos are represented, for a 
fixed value of m and a and for different values of the central density, in Fig. [3j respectively. 
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FIG. 3: General relativistic density profiles (in units of 10 25 g/cm 3 ) (left figure) and mass profiles (in units of 1O 1O M0) (right 
figure) for a Bose-Einstein condensate galactic halo with m — 5.6 x 10 -34 g, a — 10 6 fm and different values of the central 
density: p$ M = 10" 25 g/cm 3 (solid curve), p$ M = 2 x 10~ 25 g/cm 3 (dotted curve), p$ M = 3 x 10~ 25 g/cm 3 (short dashed 
curve) and p^ M = 4 x 10 -25 g/cm 3 (long dashed curve). 

The tangential velocity of a test particle moving in the dark halo is given by 



Vtg = 



\ 



v 1- 



2 M<: 



(64) 



and is represented as a function of the distance to the galactic center r in Fig. [4] 

For a galactic dark matter halo with mass of the order of Mdm — 1O 12 M0 and radius Rdm = 100 kpc, the quantity 
2GMdm/c 2 Rdm is of the order of 9.6 x 10~ 7 , which is much smaller than one. Therefore in the nominator of Eq. 
([50)1 one can neglect the term 2GmDM{r)/c 2 r. For central densities of the galactic halos of the order of 10 -26 g/cm 3 

and for a = 10 6 fm, the quantity Ao = Uop^ M /c 2 w 10~ 7 , which implies \q9 <C 1. 

Hence for galactic halos consisting of Bose-Einstein condensates, we can obtain a generalized Lane-Emden equation 
that contains the 1/c 2 corrections 



1 d 



r) 2 drj \ drj 



,d0 



= A £ 



drj ^ 9 i^drj 



(65) 



In the zeroth order of approximation (Ao — > which is equivalent to c 2 — » 00) the resulting Lane-Emden equation 
reduces to the Newtonian limit Eq. ([23]) with n = 1, whose solution is given by Eq. ([29]) . Like in the Newtonian case 
one can obtain a series expansion of the generalized Lane-Emden equation. 



VI. COMPARING THE MODEL WITH THE OBSERVATIONAL DATA 



In order to test our results we compare the predictions of our model with the observational data on the galactic 
rotation curves, obtained for a sample of low surface luminosity and dwarf galaxies in [39[ and [ic}, respectively. 
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FIG. 4: Tangential velocity of a test particle in a general relativistic Bose- Einstein condensate galactic halo with m = 5.6 x 10 -34 
g, a = 10 G fm for different values of the central density: P^m ~ 1(T 25 g/cm 3 (solid curve), p^ M = 2 x 1CT 25 g/cm 3 (dotted 
curve), p^ M = 4 x 10 _2 ° g/cm 3 (short dashed curve) and p^M — 6 x 10 -25 g/cm 3 (long dashed curve). 



Generally, in a realistic situation, a galaxy consists of a distribution baryonic (normal) matter, consisting of stars 
of mass M stari ionized gas of mass M gas , neutral hydrogen of mass Mm etc., and the dark matter of mass MnMi 
which we assume to be in the form of a Bose-Einstein condensate. The total mass of the galaxy is therefore M tot = 
M star + M gas + M m + M DM + ■■■ = M%* + M DM , where Mjf = M star + M gas + Mm + ... is the total baryonic 
mass in the galaxy. 

As for the baryonic mass, we assume it is concentrated into an inner core of radius r c , and that its mass profile 
Mfl(r) can be described by the simple relation 

/ \ 3/3 

M B (r) = Mjf — — , (66) 
\r + r c J 

where (3 = 1 for high surface brightness galaxies (HSB) and (3 = 2 for low surface brightness (LSB) and dwarf galaxies, 
respectively [f| • Therefore the tangential velocity of a test particle in our model, in which dark matter is in a form of 
a Bosc-Einstcin condensate, is given by 



v 2 tg (km/s) = 4.45 x 10' 



M*f 4 (lO lo M ) ( r \ 3P 
r(kpc) 



+ 8.06 x ptglf (l0" 25 g/cm 3 ) x (R kpc) 2 x 



sin(7rr/i?) /wr\ 
cos — 



irr/R \RJ 



(67) 



In order to test the validity of Eq. (|67|). which represents the basic prediction of our model, we have performed 12 
galaxy rotation curves fits to a sample of low surface luminosity (LSB) galaxies and dwarf galaxies, respectively. The 
observational data have been taken from [3^| for the LSB galaxies and from [4(| for the dwarf galaxies. The fitted 
rotation curves are represented by a solid curve in Fig. fSJ respectively, where the points represent the observational 
data. In all the considered cases we have adopted for (3 the value (3 = 2. The numerical results and the obtained 
values of the fitting parameters are summarized in Table fl] 

As one can see from the figure and from Table U for 9 over 12 cases there is an overall very good agreement 
between the data and the best fit curve. This strongly suggest that our model may be relevant for obtaining a correct 
description of the dark matter and its properties. A purely Newtonian description of these observational data by 
assuming that the entire matter is in the form of the baryonic matter concentrated in a core region is impossible. The 
Newtonian potential of the core is asymptotically decreasing, but the corrected rotation curve is much higher than 
the Newtonian one, thus offering the possibility to fit the rotation curves. On the other hand since LSB galaxies are 
usually considered to be dark matter dominated, reproducing their rotation curves represents a significant evidence 
in favor of the model. Moreover, fitting to rotation curves allows to determine the theory parameters, and especially 

(c) 

R, the radius of the dark matter halo, and the central density of the dark matter, p DM - The values obtained from 
the fitting are physically reasonable, indicating a radius of the galactic halos in the range of 8 — 25 kpc. On the other 
hand, the values of the central density of the dark matter (at the galactic center) predicted by our model are of the 
order of 10 -25 g/cm 3 , a value much smaller than the value of pdm (100 pc) w 10M Q /pc 3 w 7 x 10~ 22 g/cm 3 estimated 
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Galaxy 


M t B ot /W 10 M Q 


R/kpc 


pL C L/(10- 25 g/cm 3 ) 


x 2 


DDO 189 


0.39 


7.88 


5.51 


0.33 


NGC 100 


0.44 


8.38 


10.39 


1.64 


NGC 1560 


0.35 


9.86 


5.82 


5.42 


NGC 4395 


0.71 


10.22 


4.81 


11.07 


NGC 5023 


0.70 


5.67 


15.13 


2.31 


UGC 711 


0.093 


14.49 


4.56 


2.76 


UGC 3371 


0.17 


11.22 


6.12 


1.92 


UGC 5005 


0.28 


24.65 


2.056 


4.24 


UGC/AGC 8067* 


0.70 


686.83 


69.25 


58.32 


UGC/AGC 9888* 


1.49 


557.39 


85.49 


36.014 


UGC 10310 


0.086 


8.47 


9.03 


0.77 


UGC/AGC 241294* 


0.94 


950.31 


54.81 


61.015 



TABLE I: Numerical values of the fitting parameters for the fitted rotation curves. Note that \ 2 is relatively small which 
suggests that our model fits the data quite well. In view of the values of \ 2 °f f fle l° w luminosity dwarf galaxies, marked by 
*, we find that the model fits the high luminosity galaxies much better than the low luminosity ones. This is rather similar to 
the standard dark matter models, which also work best for high luminosity galaxies. 



If our model can fit quite well the observational results for the LSB galaxies, the fits for the dwarf galaxies based on 
the data from [40( require an unrealistically large extension of the dark matter, and they are statistically not significant. 
A possible explanation may be related to the peculiar behavior of the observational data near the origin: the tangential 
velocities do not tend to zero for r — ► 0, but are zero for a value of r of the order of 0.5 kpc, v tg (0.5 kpc) w 0. Thus, 
a possibility for a better fitting of the observational data for these galaxies would be to re-scale the radial coordinate 
by introducing a new parameter ro so that r — > r — ro and v tg (ro) = 0. However, the necessity of introducing a new 
parameter seems not to be physically motivated. 

Conventionally, galactic dark matter is modelled by an isothermal sphere in hydrostatic equilibrium, having two 
free parameters. These parameters are adjusted for every galaxy individually. One severe problem in that approach 
is that it seems that the rotational velocity of the last data point scales like the luminosity. Also, the Tully-Fisher 
relation L ~ v^ ut , a » 4, where L is the luminosity (in units of 1O 1O L0) and v ou t is the velocity at the outermost 
observed radial position, which holds for the high luminosity spiral galaxies [l[, actually suggests that the luminous 
matter and the related Newtonian effects must play an important role in determining the rotational velocities. 



VII. LIGHT DEFLECTION AND LENSING BY BOSE- EINSTEIN CONDENSATE DARK HALOS 



One of the ways we could in principle test the galactic dark matter model obtained in the previous Section would be 
by studying the light deflection by galaxies, and in particular by studying the deflection of photons passing through 
the region where the rotation curves are flat. Let us consider a photon approaching a galaxy from far distances. We 
will compute the deflection by assuming that the metric is given by Eq. Q470. together with Eqs. (|60p - (p^)) . 

The bending of light by the galactic gravitational field results in a deflection angle A(f> given by 



(A^) ss = 2|^(ro)-^ 00 |-7r ! 



(68) 



where (f>oo is the incident direction and ro is the coordinate radius of the closest approach to the center of the galaxy. 
Generally [36| 



4> (ro) - 4>c 



MO 



ro 



1 -1/2 



- 1 



dr 
r 



(69) 



In the dimensionless variables introduced in Eqs. (|56[) we have 

2M (r,)l- 1/2 



</> (ro) ~ </>c 



I/O 



1 - 



V 



-1/2 



(70) 
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FIG. 5: Parametric fits of a sample of 12 galactic rotation curves from the data presented in [3S| and [4(j • The points represent 
the observational values. The solid curve is the rotation curve determined with the use of Eq. (|67[1 ■ The dotted curve represent 
the tangential velocity of the Bose-Einstein condensate, while the dashed curve is the Newtonian galactic rotation curve, 
corresponding to the presence of the normal baryonic matter in the galactic core. 



We shall compare the values of the deflection angle A with a semi-realistic model for dark matter, in which it is 
assumed that the galaxy (the baryonic matter) is embedded into an isothermal mass distribution (the dark matter), 
with the density varying as p = er^/27rGr 2 , where a v is the line of sight velocity dispersion [l[. In this model it is 
assumed that the mass distribution of the dark matter is spherically symmetric. In fact, if the rotation curve is flat, 
then the mass distribution must be that of the isothermal sphere, for which we also have v tg = \/2a v . The surface 
density £ of the isothermal sphere is X (r) = a?j2Gr. For this dark matter distribution the bending angle of light is 
constant and is given by {A<f>) DM — 2-Kvf g /c 2 [42|. 

To compare the values of the deflection angle in the Bose-Einstein condensate and isothermal sphere dark matter 
models we introduce a dimensionless parameter A defined as 

(A0) DM 2ir {v tg /cf 

The variation of A as a function of the impact parameter tq is represented, for different values of the central density 
of the Bose-Einstein condensate dark matter, in Fig. [5] 
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FIG. 6: The parameter A as a function of the impact parameter rja = ra/r* for a Bose-Einstein condensate with m = 5.6 x 10~ 34 
g and a — 10 6 fm, for different values of the central density: p^ M = 10" 27 g/cm 3 (solid curve), p^ M = 2 x 10 -27 g/cm 3 (dotted 
curve), p$ M = 4 x 10~ 27 g/cm 3 (short dashed curve), p^ M = 6 x 10~ 27 g/cm 3 (dashed curve), p^ M = 8 x 10~ 27 g/cm 3 (long 
dashed curve) and p^ M = 10~ 26 g/cm 3 (ultra-long dashed curve). In each case vt g = 300 km/s. 




Once the light deflection angle is known, one can study the gravitational lensing by the dark matter halos. The 
lensing geometry is illustrated in Fig. [7l 

The light emitted by the source S is deflected by the lens L (a galaxy in our case) and reaches the observer O at 
an angle 9 to the optic axis OL, instead of j3. The lens L is located at a distance Dl to the observer and a distance 
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FIG. 7: The lensing geometry, showing the location of the observer O, of the lensing galaxy L and of the source S. The 
deflection angle is Acf). The angular diameter distances Dl, Dls and Ds are also indicated. 



Dls to the source, respectively, while the observer-source distance is Ds- ro is the impact factor (distance of closest 
approach) of the photon beam. 
The lens equation is given by [Hf 



tan /3 = tan0 - [tan 9 + tan ((A</)) BE -9)}. (72) 
Ds 

By assuming that the angle 9 is small, we have tan 9 « 9 and the lens equation can be written as 



p K e-^{A4>) BB - (73) 

L>S 



In the special case of the perfect alignment of the source, lens and observers, P = 0, and the azimuthal axial 
symmetry of the problem yields a ring image, the Einstein ring, with angular radius 

9 (BE)^DLS m (?4) 

J->s 

This equation can be expressed in a more familiar form by taking into account that the impact parameter ro ~ D^9, 
which gives 



where 9^^ is the angular radius of the Einstein ring in the case of standard general relativity, 9^^ = 
\pk {Dls/DsDl) GMb, with Mb denoting the baryonic mass of the galaxy. 

In the case of a galaxy with a heavy isothermal dark matter distribution, the Einstein radius of the lens formed in 
perfect alignment is [13] 

eT ] - (^) (76) 



18 



The ratio A of the Einstein's rings angular diameters in the brane world models and in the isothermal dark galactic 
halo model is 

Therefore the ratio of the angular radii of the Einstein rings in the Bose-Einstein condensate dark matter model 
and in the isothermal dark matter model is given by the same parameter A which has been already introduced in Eq. 
(I7T1) . Hence the variation of the ratio of the Einstein rings in the two models can also be obtained from Fig. [5J 



VIII. DISCUSSIONS AND FINAL REMARKS 



Galactic rotation curves pose a challenge to present day physics and one would like to have a better understanding 
of some of their intriguing phenomena, like their universality and the very good correlation between the amount of 
dark matter and the luminous matter in the galaxy. In the present work we have considered, and further developed, 
an alternative view to the dark matter problem, namely, that the galactic rotation curves can be explained by models 
in which dark matter is in the form of a Bose-Einstein condensate. This assumption leads to a complete description of 
the basic properties of the dark matter condensates at both the Newtonian and general relativistic level, and gives a 
definite set of predictions which can be tested observationally. In particular, we have shown that the model provides 
a good descriptions for the galactic rotation curves by fitting our model to the data of 12 observed curves. 

The numerical values of the basic parameters (mass and radius) of the condensed dark matter halos sensitively 

(c) 

depend on the mass m of the condensate, on the scattering length a and on the central density p DM , Rdm = 

R{a, m, p*d M ), Mjjm = Mdm(o>, m, Pdm)- Of course, in general, the values of the mass and radius of the gravitational 
condensate depend on the adopted model for the non-linearity. 

The scattering length a is defined as the zero-energy limit of the scattering amplitude / [12j . Depending on the 
spin dependence of the underlying particle interaction, the scattering length may in general be also spin dependent. 
The spin independent part of the quantity is referred to as the coherent scattering length a c . The scattering lengths 
can be obtained for some systems in the laboratory. In our estimations we have used the value of the scattering 
length obtained for dilute atomic Bose-Einstein condensates in terrestrial laboratory experiments. Another essential 
parameter is the mass m of the condensate particle, which, due to the lack of information about the physical nature 
of the dark matter, is a free parameter, which must be constrained by observations. 

A powerful observational tool for discriminating between standard dark matter and Bose-Einstein condensate models 
is the study of the deflection of light (gravitational lensing) by the dark matter halos. Due to the fixed form of the 
galactic metric in the flat rotation curves region, in standard dark matter models the light bending angle is a function 
of the tangential velocity of particles in stable circular orbit and of the baryonic mass and radius of the galaxy. 
Generally, the specific form of the bending angle is determined by the galactic metric, and this form is very different 
for the Bose-Einstein condensate dark matter as compared to the other dark matter or modified gravity models 
(MOND 3, non-symmetric gravity ||, long-range self-interacting scalar fields [ll[ etc.). 

The gravitational light deflection angle predicted by the Bose-Einstein condensate models is much larger than the 
value predicted by the standard dark matter approach. For example, there are significant differences in the lensing 
effect with respect to the isothermal dark matter halo model. Therefore the study of the gravitational lensing may 
provide evidence for the existence of the dark matter in the form of a Bose-Einstein condensate. Therefore, the 
study of the galaxy-galaxy lensing and of the dark matter halos' properties could provide strong constraints on the 
Bose-Einstein dark matter model and on related physical models. 

In the present approach to dark matter all the relevant physical quantities can be obtained from observable pa- 
rameters (the dark halo mass, the radius of the galaxy and the observed flat rotational velocity curves). Therefore 
this opens the possibility of testing the Bose-Einstein condensation models by using astronomical and astrophysical 
observations at the galactic-intergalactic scale. In this paper we have provided some basic theoretical tools necessary 
for the in depth comparison of the predictions of the condensate model and of the observational results. 
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APPENDIX A: MOTION OF TEST PARTICLES 

The motion of a test particle with four- velocity in the gravitational field of the galaxy can be described by the 
Lagrangian 0] 

where d£l 2 = d0 2 + sin 2 9d<p 2 and dr — cdt. From the Lagrange equations it follows that we have two constants 
of motion, the energy E = e u ( r 'i and the angular momentum I = r 2 4> 37]. The condition u^u^ = —1 gives —1 = 
—e"( r H 2 + e A ( r V 2 + r 2 (f> 2 and, with the use of the constants of motion we obtain 

E 2 = e » +x r 2 +e"(J? s + lj. (A2) 

This equation shows that the radial motion of the particles on a geodesic is the same as that of a particle with 
position dependent mass and with energy E 2 /2 in ordinary Newtonian mechanics moving in the effective potential 
V eff (r) = e v{r) (l 2 /r 2 + l). The conditions for circular orbits dV eff /dr = and r = lead to [S3] 



1 rV 



r = - -,E* = r- (A3) 

2 ^ Tv \ rv 

The line element given by Eq. (|47p can be rewritten in terms of the spatial components of the velocity, normalized 
with the speed of light, measured by an inertial observer far from the source, as ds 2 — ~c 2 dt 2 (l — v 2 /c 2 ) 7], where 



For a stable circular orbit r = 0, and the tangential velocity of the test particle can be expressed as 

c 2 e v \cdt) 

In terms of the conserved quantities the angular velocity is given by 

c r 2 E 2 

With the use of Eqs. (|A3|) we obtain 



(A4) 



(A5) 



(A6) 



